oo 



> 



o\ 



X 



ITFA-97-02 
February 1997 



Numerical study of plasmon properties in the 

SU(2)-Higgs model 
S' 

^ '. Wai Hung Tang* and Jan Smit^ 

bJQ! Institute of Theoretical Physics, University of Amsterdam 

<£ ■ Valckenierstraat 65, 1018 XE Amsterdam, the Netherlands 



revised version, July 1997 



Abstract 



o 

r—N | Using the (effective) classical approximation, we compute numerically 

t^- ■ time-dependent correlation functions in the SU(2)-Higgs model around the 

electroweak phase transition, for m# ss mw- The parameters of the clas- 
sical model have been determined previously by the dimensional reduction 
relations for time-independent correlators. The H and W correlation func- 
tions correspond to gauge invariant fields. They show damped oscillatory 
behavior from which we extract frequencies to and damping rates 7. In 
the Higgs phase the damping rates have roughly the values obtained in 
analytic calculations in the quantum theory. In the plasma phase (where 
analytic estimates for gauge invariant fields are not available), the damp- 
ing rate associated with H is an order of magnitude larger than in the 
Higgs phase, while the W correlator appears to be overdamped, with a 
small rate. The frequency ujh shows a clear dip at the transition. The 
results are approximately independent of the lattice spacing, but this ap- 
pears to be compatible with the lattice spacing dependence expected from 
perturbation theory. 
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1 Introduction 

The classical approximation has been introduced some time ago in quantum field 
theory to avoid the difficulties of numerical simulations of time dependent phe- 
nomena at high temperature [|IJ. It is potentially a very powerful method for 
nonequilibrium phenomena as well. It has been used for the numerical study of 
the Chern-Simons diffusion rate in the (1+1)D abelian Higgs model [2-7], the 
(3+l)D SU(2)-Higgs model H i, and pure SU(2) gauge theory 0, as well 



as in a study of real time properties of the electroweak phase transition [12]. The 
results have been encouraging. For example, the sphaleron rate in the classical 
abelian Higgs model obtained by numerical simulations agrees with the analytic 
expression of quantum field theory [4-7]. In the corresponding pure SU(2) and 
SU(2)-Higgs cases the rate appeared to be lattice spacing independent in refs. 
| [PI] , |10|1 , but these results are now considered to be misleading, because of bad 
ultraviolet properties of the observable used for estimating the Chern-Simons 
diffusion. Recent simulations with improved observables have led to much bet- 
ter results in the Higgs phase and somewhat different values for the rate in the 



plasma phase [[U| [HJ. The case of SU(3) gauge theory has also been studied ][L5[ . 
However, these new results for the rate show lattice spacing dependence. 

A recent analysis of the dynamics of hot quantum gauge fields led to the 



conclusion that the classical rate should depend on the lattice spacing [|T6j . Other 
lattice artefacts such as anisotropies have been anticipated earlier |T7[. These 
have been studied in detail in |Pj|, where a useful fudge factor was derived for 
converting the rate in the classical theory to an estimate for the rate in the 
quantum theory. These studies were analytic, based on a separation of scales 
gT, g 2 T, in resummed perturbation theory. This may be questionable when the 
gauge coupling g is not very small. 

In any case, the very validity of the effectively^] classical approximation for 
real time processes needs further clarification. For this reason we undertook the 
present study of correlation functions in the classical SU(2)-Higgs theory. The 
emphasis here is on time dependent correlation functions. The static correla- 
tion functions, which describe the initial conditions of the time dependent ones, 
have been well studied since they correspond to a dimensional reduction approx- 
imation of the quantum theory at high temperature. The static theory can be 
renormalized, with finite renormalizations determined by matching the classical 
to the quantum theory ||19| . An important question is if the time dependent 
theory is also renormalizable. 

Recently |[20|| , time dependent correlation functions have been studied in hot 
classical 4 theory using perturbation theory and it was shown that the ultraviolet 
divergencies can be absorbed in counterterms, the same mass counterterms as 



1 By 'classical' we shall always mean 'effectively classical' in this paper, i.e. in terms of a 
classical action with effective coupling constants. 



needed in the static theory. This led in particular to a finite classical plasmon 
damping rate. Matching the parameters of the classical static theory with those 
in the quantum theory, this classical damping rate was found to agree with the 
hard thermal loop expression of the quantum theory. Such matching was applied 
previously to the SU(2)-Higgs model for the computation of the Chern-Simons 
diffusion rate |10| . The analysis is to be completed by a more detailed comparison 
between the fully time dependent classical and quantum theories [pl| . 

Our approach is different from propositions based on an effective theory of 
hard thermal loops |L7], 22], in which the cutoff of the classical theory should 



not be removed. We are instead aiming at the construction of a classical SU(2)- 
Higgs theory that is renormalizable (as in the scalar case fl20|), such that physical 



results are free of regularization artefacts, for sufficiently large cutoff or small 
lattice spacing. We do not know yet if such a construction exists. It may involve 



nonlocal counterterms or other degrees of freedom, e.g. as in pB| |. The reason 
is that in gauge theories the divergences in the classical theory are non local, 
of a form that is analogous to the hard thermal loops in the quantum theory. 
Furthermore, the form of the usual classical hamiltonian does not allow for a 



counterterm to compensate for the divergence in the Debije mass [ 10 1 . Also the 



plasmon frequency is divergent fL8fl . On the other hand, in the static case, the 
large Debije mass is the rationale for integrating out the corresponding modes to 
get a simpler form of the dimensional reduction approximation. Similarly, in the 
time dependent case, without additional nonlocal counterterms, any modes with 
divergent frequencies or damping rates may simply decouple from coarse grained 
large time behavior. The question is, if what is left over can be interpreted as a 
(nonperturbatively) renormalizable theory. 

In this paper we deal with the usual classical model and concentrate on rela- 
tively simple time dependent correlation functions. First of all, we want to see if 
it is feasible to compute such autocorrelation functions numerically with reason- 
able accuracy, find out how the results depend on the lattice spacing, using the 
counterterms known from the static theory, and compare with hard thermal loop 
results of the quantum theory. 

Such a comparison is not straightforward in the high temperature phase. In 
the analytic calculations the correlation functions are constructed from the ele- 
mentary field variables using gauge fixing, while nonperturbatively one tends to 
use gauge invariant composite local fields. Even if the final results obtained ana- 
lytically are gauge invariant, it is not clear that these refer to the same quantities 
as those obtained with the gauge invariant fields. For example, in the plasma 
phase there is a large discrepancy between the two cases in the screening masses 



with W quantum numbers [24-29] (see also the reviews JJI], |3T[] and the inter- 
pretation given in ref. ||32|| ). In the Higgs phase the results agree, as can be 
understood from the usual argument which we now review. 

The simplest gauge invariant composite fields (which we also use) read in the 



continuum 




h=<pv ?=r B > w 



W£ = zIV0t( 9 _ zA%)0t a , = ^/« , (2) 

y-fufdj 

with yj the Higgs doublet, A^ the gauge fields and t a = r Q /2 the SU(2) generators. 
In the Higgs phase these fields correspond to the usual Higgs and gauge fields. 
Since 4> = VHV with V G SU(2), W^ is essentially the gauge field in the unitary 
gauge V — 1, as can be seen in perturbation theory by replacing (p by its vacuum 
expectation value, neglecting fluctuations. This is not useful in the plasma phase 
where the theory behaves like QCD and perturbation theory runs into infrared 
problems. Nonperturbatively H and W^ (m = 1, 2, 3) are primarily characterized 
by their behavior under spatial rotations and weak isospin: (scalar, scalar) and 
(vector, vector), respectively. We have not studied the time component Wq, which 
transforms as (scalar,vector). 

To get an idea of the properties of the system being simulated we also com- 
puted static correlators of H and W and estimated their screening masses. The 
temperatures varied around the electroweak transition and parameters were fur- 
thermore chosen such that for the zero temperature masses mg ~ mw From 
the time dependent correlation functions we then attempted to extract plasmon 
frequencies and damping rates of H and W. 

We should warn the reader about our terminology. For convenience we often 
call in this paper the frequencies and attenuation rates associated with the com- 
posite W and H fields at zero momentum: 'plasmon' frequencies and damping 
rates. This is controversial, since the terminology 'plasmon' is normally associ- 
ated with the collective modes of the gauge field A^. 

2 Classical SU(2)-Higgs model at finite temper- 
ature 

We use the same notation as in [|K|]. The classical SU(2)-Higgs model is defined 
on a spatial lattice with lattice distance a. The parallel transporters (lattice 
gauge field) are denoted by t/ mx , D m (p^ = t/ mx (^ x+ ^ l — <^ x is the covariant lattice 
derivative acting on the lattice Higgs doublet </5 x , U mnil is the product of par- 
allel transporters around the plaquette (xmn). The canonical momenta in the 
temporal gauge are denoted by 7f, 7ft and E^, with nontrivial Poisson brackets 

{U mx , E^ x } = iU mx r a , {^ x , vr x } = 1. (3) 



The effective classical hamiltonian is given by 

^ = $ B > ( 4 ) 

x m<n 

+ W^x + Afe-D 2 ) 2 ]. (5) 



where T is the temperature and (3, A and v 2 are effective couplings which depend 
on a, T and the couplings in the 4D quantum theory, 



4 - m 



/?--^ A«-S-, (6) 
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4Au 2 = a 2 T 2 



T 2 V2 ) V 18 3 aT ' ' ' 

3 \ 4 1 /149 + 9p 27 + 6p - p 2 , aT 27r? + 6pr7 - p 2 r? 
„J n — ^7 — + ™ ln ^r 



2*/ 8tt 2 V 486 27 2 27 

£ = 0.252731. (8) 

Here p = mf^/m 2 ^ and g and mj are parameters in the MS-bar scheme at scale 
fiT = 47rexp(— 7b)T ~ 7T. These relations follow from matching the static 
classical theory to the dimensionally reduced quantum theory [f[T[j , using the 
results of ref. [Q. As in JL0|| we set ze = z n = I. The implementation of these 



parameters has to wait for a detailed matching between the time dependent 
classical and quantum theories. 

The equations of motion follow from the above hamiltonian and Poisson brack- 
ets. The initial conditions are distributed according to the classical partition 
function 

Z = J DEDUDnDcp [J] 6(G«)} exp(-# cff /T), (9) 

xct 

where S(G^) enforces the Gauss constraint, which is part of the temporal gauge 
formalism. Time dependent correlation functions of observables O are defined by 
averaging O(t)O(0) over initial configurations, 

(O(t)O(O)) = ^fDEDUDirD<p[Y[6(GZ)] exp(-# eff /T) O(t)O(0). (10) 

xa 

Examples of such observables are the simple gauge invariant fields we use to study 
the W and Higgs excitations: 

WZ x =iTr<&U m3t (l> x+A ,Ta, (11) 

H x = yfe- (12) 



ft 


N 


Ph 


vl 


T c /m H 


12 


20 


0.347733 


0.26295 


2.14 


12 


24 


0.34772 


0.26273 


2.15 


20 


32 


0.34173 


0.15597 


2.15 



Table 1: Critical values of (ft,v) for A = 1/2. The values of T c /rriH were 
calculated using eq. (0). 



The static system (i.e. the system described by the classical partition function 
at time zero) is believed to be well understood because of its connection with 
dimensional reduction. For given v 2 and relatively small A there is a critical 
value of (3, above which the system is in the Higgs phase and below which it is 
in the plasma phase. The transition is of first order |34| , weakening in strength 
as m H /mw increases, and turning into a crossover near m H « mw p5 |. 

As in ||10|| , we have chosen A = 1/2, which means vnnjra^ ~ 1. The transition 
is then still clearly visible, especially in the rate of Chern-Simons diffusion [[KJ . 
Accordingly set p = 1 in the application of eq. (0). We also neglect the running 
of g and mj with temperature, and choose g = 2/3, which implies aT = 9/J3. 
For given v 2 eq. (|?p then gives T/mn as a function of /3. 

Keeping instead T/tuh fixed while increasing j3 and changing v 2 according to 
(0) should get us closer to the continuum limit. Physical quantities of the static 
theory then become [3 independent. It is not clear at this point if the same holds 
for physical quantities related to the full time dependent theory, such as plasmon 
masses and rates. 

Table [l] shows some pseudo critical values (3 C and v 2 on lattices of size N 3 . 
These values are obtained by a translation of the dimensional reduction results 
of [34| into the parameterization used here in terms of A and v 2 . For the record 
we mention the connection with the parameterization used in [34"|: /3 = j3c, 
A = WoPr/Ph, (^VA) = 4(3/?# - 2{3 R - 1)/(3 H , ^ = ^ = 0. The last two 
parameters were not zero in the action used in [5^] , but the difference does not 
seem to be important for the location of the phase transition ||10| . 



3 Numerical computation 



We use the algorithm and numerical implementation offered in ref. [33]. A simula 



tion consisted of alternating Langevin runs (the generation of initial conditions) 
and Hamilton runs (integration of the equations of motion). Spatial correla- 
tion functions used for the study of screening properties were averaged over the 
Langevin runs. The Langevin runs lasted typically 120 (in lattice units), which 
we found sufficient to decorrelate observables such as the 'link' Ti<f)^U<j), while 
the Hamilton runs lasted in most cases several thousand. Details on the statistics 



will be given below where we summarize the parameter values in Table 0. 

The spatial correlation in the 3-direction at distance d, of a local observable 
O x at zero transverse momentum, was estimated as follows, 

°o( d ) = jj E [Ox 3 +dO X3 - d~^ d oZ\ , (is) 

£3 
7T7^ / , ^Xi,X2,X3; (,-1-4) 



^£3 



/AT2 

£1,2:2 



where the bar denotes the average over the Langevin bins (we used periodic spatial 
boundary conditions). In the limit of an infinite number of bins O x — > (O x ) 
which is independent of x. For the Higgs mode we used O x = H x (cf. (|12j) 
and we averaged the correlators over the three different directions (3, 1, and 2). 
For the jy-mode we used the transverse components O x = J2 a =i ^ibo k — 1,2, 
which leads to a Gw k (d) independent of k, and we averaged over k, Gw(d) = 
\Gw\{d) + Gw 2 (d)]/2, and over the three different directions. (The summation 
over a = 1, 2, 3 picks out a particular direction in isospin space; we could have 
improved statistics a little by correlating and averaging the a's as well). 

For the time dependent correlators we used microcanonical averaging over the 
Hamilton parts of the simulation in addition to canonical averaging over initial 
conditions at the end of the Langevin parts. The autocorrelation functions were 
constructed analogously to (|13"1), 

1 f tb 



Co{t) = - / dt Q O(t +t)O(t )-O(t + t) O(t ) , (15) 

h Jo l J 

where < t < t run — £&. Here the bar denotes the canonical average and t run is 
the maximum time in the Hamilton run. When the number of initial conditions 
goes to infinity this expression approaches the exact (O(t)O(O)) — (O) 2 . For 
the determination of the plasmon properties we used for O the zero momentum 
projections 

V XQ V X 

and for W we averaged over k. 



4 Results for the screening correlators 

Table |2] gives a summary of the parameter values used in our simulations. The 
values of T/tuh correspond to eq. (0). (For f3 — 6, 8 the lattice spacing is so 
large that eq. (|7|) breaks down; note that 0(a) terms are neglected in [■ • ■].) The 
ratios T/T c have been added for convenience and follow simply by devision by 
T c / Mh(ht c ) in Table [I], i.e. without taking into account the running of mainr) 



p 


N 


v 2 


T/m H 


T/T c 


t-run 


nc(sc) 


nc(pl) 


20 


24 


0.263 


0.88 


0.41 


1000 


15 


15 


14 


24 


0.263 


1.59 


0.74 


6400 


15 


18 


21.7 


32 


0.156 


1.67 


0.78 


7000 


31 


12 


13 


24 


0.263 


1.82 


0.85 


5000 


42 


42 


13' 


20 


0.263 


1.81 


0.85 


5000 


40 


11 


12.5 


24 


0.263 


1.98 


0.92 


6400 


35 


30 


12 


24 


0.263 


2.15 


1 


1500 


30 


29 


11 


24 


0.263 


2.64 


1.23 


400 


89 


86 


11' 


20 


0.263 


2.63 


1.23 


2500 


57 


56 


11* 


20 


0.246 


3.26 


1.52 


3000 


61 


61 


18.3 


32 


0.156 


3.26 


1.52 


2500 


37 


21 


10 


24 


0.263 


3.46 


1.61 


1280 


24 


48 


8 


24 


0.263 


- 


- 


320 


41 


51 


6 


24 


0.263 


- 


- 


640 


35 


22 



Table 2: Summary of parameter values and statistics, ordered according to 
increasing T/T c . The prime or * on j3 distinguishes different N and v 2 . 



with T. The time t run is the length of the Hamilton runs used in the computation 
of the autocorrelation functions, in lattice units; nc(sc) and nc(pl) are the number 
configurations (number of Langevin runs) used in the calculation of the screening 
and plasmon properties, respectively. 

We first made a scanning simulation on the N = 24 lattice with varying 
(3 and v 2 fixed at 0.263, its value for /3 C = 12. This means the temperature 
varied as shown in Table 0. By fitting the correlation functions GH,w{d) to 
Acosh[am(d — N/2)] we obtained the screening masses m shown in Table |[ 

We made a least squares fit to the data at d — 3 - 5, which led to reasonable 
looking fits, and the errors were obtained with the jackknife method. This pro- 
cedure is of course not state of the art as in |28[ , but we only wanted to get an 
overall impression of the data at the chosen parameters using a moderate amount 
of computational effort. 

After this scan we decided to concentrate on two values of /3 corresponding to 
two values of the lattice spacing, in the Higgs phase as well as in the plasma phase. 
The parameters v 2 and N were changed accordingly such that the two (3 values 
tentatively describe the same physical situation. The guiding parameter values 
were those of the transition obtained nonperturbatively in [3~3], which are shown in 
Table||: (/?, N, v 2 ) = (12, 20, 0.263) and (20, 32, 0.156). Modulo scaling violations 
this corresponds to a ratio of lattice spacings 12/20 = 0.6 at approximately the 
same physical volume (20/32 = 0.625). The new parameters in the Higgs phase 
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p 


iV 


T/T c 


arriH 


amw 


m H /g 2 T 


m w /g 2 T 


20 


24 


0.41 


0.59(5) 


0.50(2) 


2.95(25) 


2.5(1) 


14 


24 


0.74 


0.34(1) 


0.41(2) 


1.19(4) 


1.44(7) 


21.7 


32 


0.78 


0.22(2) 


0.25(2) 


1.19(11) 


1.36(11) 


13 


24 


0.85 


0.28(3) 


0.35(2) 


0.91(9) 


1.14(7) 


13' 


20 


0.85 


0.26(2) 


0.39(3) 


0.85(7) 


1.27(10) 


12.5 


24 


0.92 


0.22(3) 


0.38(1) 


0.69(9) 


1.19(3) 


12 


24 


1 


0.146(29) 


0.37(2) 


0.44(9) 


1.11(6) 


11 


24 


1.23 


0.46(3) 


0.88(14) 


1.27(8) 


2.4(4) 


11' 


20 


1.23 


0.44(3) 


0.96(11) 


1.21(8) 


2.64(30) 


11* 


20 


1.52 


0.59(5) 


1.06(20) 


1.62(14) 


2.92(55) 


18.3 


32 


1.52 


0.31(3) 


0.65(10) 


1.42(14) 


2.97(46) 


10 


24 


1.61 


0.68(4) 


1.34(28) 


1.7(1) 


3.4(7) 


8 


24 


- 


1.12(20) 


1.04(58) 


2.24(40) 


2.1(12) 


6 


24 


- 


1.2(6) 


1.5(18) 


1.8(9) 


2.3(27) 



Table 3: Results for the screening masses. 



were obtained by changing (3 = 12 — > 13 and ft = 20 — > 21.7 ~ 13/0.6, keeping 
fixed the corresponding N and v 2 . These pairs of (j3,N,v 2 ) values correspond 
roughly to the same T/M^ (using the matching formula (]?]), as shown in Table 
(cf. (3 = 13' and 21.7). The same parameter combinations were used in JT(| for 
the computation of the Chern-Simons diffusion rate. Similarly the plasma phase 
parameters were obtained by changing (3 — 12 — > 11 and /3 = 20 — > 18.3 ~ 11/0.6, 
again without changing the corresponding N and v 2 . However, here the matching 
formula (|7]) gives a larger mismatch in T/rnn (cf. /3 = 11', 18.3 in Table ^j. To 
keep T/ttih fixed more precisely we have to adjust v 2 . Therefore we also did 
measurements in the plasma phase at ((3,N,v 2 ) = (11,20,0.246), for which v 2 
was chosen such that T/m H equals its value at (j3,N,v 2 ) = (18.3,32,0.156). In 
the Tables the J3 value is denoted by 11*. 

Table |3| also shows the results for the screening masses of these additional 
simulations. We fitted the data at d = 3 - 5 for iV = 20 and at d = 4 - 7 
for iV = 32. Fig. |l] shows the screening masses in physical units, m/g 2 T, versus 
T/T c . The transition at T = T c is clearly visible and we also see the ratio m^/m^ 
approaching 1 at low temperatures, as expected from the choice A = 1/2. (The 
dimensional reduction approximation is of course questionable at the most left 
points where T/Mh is only 0.88.) We observe no significant volume dependence 
away from the phase transition in the data at T/T c = 0.85 and 1.23. The data 
at T/T c = 1.52 (20 3 , (3 = 11* and 32 3 , 18.3) show reasonable lattice spacing 
independence (recall that the ratio of lattice spacings is 0.6). The difference in 



m/g 2 T between (3 = 13' (20 3 ) and 21.7 (32 3 ) is evidently due to the difference in 
physical temperatures, T/T c = 0.84 and 0.78, and not a signal of lattice spacing 
dependence. 



p 


am e sH 


am c gw 


m eSH /g 2 T 


m eSW /g 2 T 


13' 


0.24(2) 


0.43(4) 


0.78(7) 


1.40(13) 


21.7 


0.18(3) 


0.27(3) 


0.98(16) 


1.46(16) 


11' 


0.64(5) 


1.65(14) 


1.76(14) 


4.5(4) 


11* 


0.64(4) 


1.84(13) 


1.76(11) 


5.1(4) 


18.3 


0.32(5) 


1.26(11) 


1.46(23) 


5.8(5) 



Table 4: Effective screening masses from eq. (JT8| 



We also like to mention our experience with a momentum space analysis of the 
correlation functions, which has been popular in numerical studies of the Higgs- 
Yukawa sector of the Standard Model. This analysis is based on the Fourier 
transform of GH,w(d), 



G H Mp) = J2 e ~ tPdG ">w(d). 



(17) 



Here p is one of the lattice momenta, ap = 27m/ N, n — 0, 1,- ■ ■ ,N — 1. Within 
statistical errors these are real functions of p. We deduced effective screening 
masses m c g from the expansion 

G{p)- 1 = Z-^mls + p 2 + 0{p% (18) 

by fitting G~ l as a function of p 2 = 2 — 2 cos ap to a straight line, using two of the 
lower p values (n = ap/(2ir/N) = 0, 1 and 1,2). The resulting effective screening 
masses (cf. Table |j) tend to be significantly larger than the true screening masses 
above, especially for W in the plasma phase, for which G~ 1 (p) shows strong cur- 
vature. (Fitting more data points to a quadratic function had no significant effect 
on the results.) The effective masses are expected to be a good approximation 
to the screening masses only in case of pole dominance of G(p), which evidently 
is not the case for W in the plasma phase. 

It is known that in the Higgs phase the simple fields ( O , |T^ ) which we use 
perform reasonably well in creating dominantly the H and W out of the vacuum, 
but in the plasma phase they easily create also excited states with the same 
quantum numbers as H and W (for a recent study see p8"|). The excited state 
contribution weakens of course pole dominance. To do better one has to use more 
sophisticated observables instead of the simple if x and Wj^, e.g. the smeared 
fields used in [J2"8 |. 
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5 Plasmon correlators in the Higgs phase 

The results in the previous section for the screening properties of H and W give 
us a feeling how well we are doing in a familiar situation, in the theoretical frame- 
work of dimensional reduction. We now turn to the time dependent correlation 
functions, for which we are on less established ground. Figs. |] and [3] show an 
initial scan of the autocorrelation functions C#(£) and Cw{t) for various j3 at 
fixed v 2 = v 2 (/3 = 12) on the 24 3 lattice. The horizontal time scale is in lattice 
units (i.e. (number of time integration steps) x (step size) as it appears in the 
computer code). The time in lattice units can be converted to physical units by 
multiplication with e.g. the temperature in lattice units. We see the period of 
the oscillations increasing as we approach the transition at ft — 12 from above 
(from the Higgs side)0. 

In the H case (Fig. U) there is a drastic change of behavior below (3 = 12 
(the plasma side), where the signal is much smaller. As the insert shows the data 
here is very noisy, although we think the very short time (t < 2) data is still 
relevant, because short times allow for more microcanonical averaging. Notice 
the oscillations in this very short time region. 

For W (Fig. ^|) the transition appears more gradual. As j3 is lowered down 
from 20, small oscillations in the small time region (t < 2) appear already in the 
Higgs phase. When j3 is decreased further below 12 the initial oscillations increase 
somewhat and so does the slope afterwards, which is nearly zero for /3 = 11, 10, 
8. In this case we have stopped the drawing of the curves at early times because 
the statistical noise would be overwhelming. 

To extract a plasmon frequency and damping rate we fit the data in the Higgs 
phase to the simple asymptotic form 

Re^ cos(cut + a) + C. (19) 

For a free scalar field we would have R = T/m 2 , 7 = C = 0. With interactions 
the form flT9| ) is expected to be a good approximation in case of small damping 



3(J. It turns out that the phase a is not really needed and may be set equal to 
zero. In the plasma phase the signal of our 24 3 data is too noisy for a quantitative 
analysis. For the Higgs phase the 24 3 data led to the plasmon parameters shown 
in Tables f| and |6|. The errors are based on jackknifing with respect to the 
initial conditions. The damping rate 7 is sensitive to the beginning of the fitting 
range, the fits above started roughly from the third maximum of Cn(t) (the first 
maximum is at t — 0). 

We increased the numerical effort on the plasma phase parameters values 
(3 = 11', 11* and 18.3 and did similar runs for 13' and 21.7. Fig. ^ shows the 
autocorrelation function C#(£) in the Higgs phase for j3 = 13' and 21.7. A fit 

2 The curves at J3 = 12 belong clearly to the Higgs phase sequence, which suggests that J3 C 
is actually slightly smaller that 12. 



11 



p 


N 


T/T c 


CLUJh 


au>w 


u H /g 2 T 


uj w /g 2 T 


20 


24 


0.41 


0.508(2) 


0.551(2) 


2.54(1) 


2.76(1) 


14 


24 


0.74 


0.350(4) 


0.430(6) 


1.23(2) 


1.51(2) 


21.7 


32 


0.78 


0.210(2) 


0.296(2) 


1.14(1) 


1.61(1) 


13 


24 


0.85 


0.28(2) 


0.38(2) 


0.91(7) 


1.24(7) 


13' 


20 


0.85 


0.284(5) 


0.382(5) 


0.92(2) 


1.24(2) 


12.5 


24 


0.92 


0.24(1) 


0.35(2) 


0.75(3) 


1.09(6) 


12 


24 


1 


0.16(1) 


— 


0.48(3) 


— 


11' 


20 


1.23 


0.33(3) 


— 


0.91(8) 


— 


11* 


20 


1.52 


0.44(3) 


— 


1.21(8) 


— 


18.3 


32 


1.52 


0.23(2) 


— 


1.05(9) 


— 



Table 5: Results for the plasmon frequencies. 



to (|19"D starting roughly at the third maximum led to the plasmon parameters 
shown in Tables |5] and |6|. Figure |5| shows the M^-autocorrelation functions in the 
Higgs phase. It is clear that the W damping rate is considerably larger than 
for H and to avoid running into noise we had to lower the starting time of the 
fit to (|T9"D down to roughly the second maximum. The constant C of the fits is 
non-zero. However, we have seen the autocorrelators to tend to zero at very large 
times. The results in physical units are shown in Figs. ^| and |7|. From the smooth 
adjustment of the (3 = 21.7 (32 3 ) data with the J3 = 13' and 13 (20 3 and 24 3 ) 
data it follows that the plasmon frequencies and damping rates are nearly lattice 
spacing independent. 

However, the 32 3 data point in Fig. ^ for the W plasmon frequency at T/T c = 
0.78 does look somewhat high compared to the trend of the 24 3 data. Interpolat- 
ing the 24 3 data linearly between T/T c = 0.74 and 0.85 gives us a hypothetical 
data point u w /g 2 T = 1.41(3) at T/T c = 0.78, J3 = 13.6. The 32 3 data point at 
T/T c = 0.78, 1.61(1), is a factor 1.14 higher. This is much less than the corre- 
sponding ratio of lattice spacings, 21.7/13.6 = 1.60, and in the previous version 
of this paper we ignored this 14% discrepancy as due to a possible systematic 
error in our data analysis. Now we shall take it seriously and try to interpret 
it as an expression of lattice spacing dependence. As a rough model, let us as- 
sume that io 2 u = 2Xv 2 , Loly = g 2 v 2 /4 + u 2 , where v 2 is simply defined by the 
equation v 2 = cj^/2A, and uj 2 is an additional contribution to the IV-plasmon 
frequency. Since the Higgs parameters are tuned according to dimensional reduc- 
tion we assume that u 2 H is independent of a. This is consistent with the trend 
of the data in Fig 
loop arguments 



]. With X/g 2 = 1/8 we get u% = ur H + uj 2 . Hard thermal 
lead to the conclusion that there is a divergent component 
in u$y, given by 0.08Q0QC A g 2 T/a, with C A = 2 for SU(2). We have checked in 
a one loop calculation perturbative calculation on the lattice that a scalar gives 
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N 


T/T c 


a'jH 


a^ w 


lH/g 2 T 


iw/g 2 T 


20 


24 


0.41 


0.0012(19) 


0.008(4) 


0.006(10) 


0.04(2) 


14 


24 


0.74 


0.0045(31) 


0.03(1) 


0.016(11) 


0.105(35) 


21.7 


32 


0.78 


0.0052(8) 


0.041(8) 


0.028(4) 


0.22(4) 


13 


24 


0.85 


0.013(29) 


0.04(5) 


0.042(94) 


0.13(16) 


13' 


20 


0.85 


0.011(5) 


0.08(2) 


0.036(16) 


0.26(7) 


12.5 


24 


0.92 


0.012(8) 


0.09(9) 


0.038(25) 


0.28(28) 


12 


24 


1 


0.03(2) 


- 


0.09(6) 


— 


11' 


20 


1.23 


0.079(61) 


0.0116(17) 


0.22(17) 


0.032(5) 


11* 


20 


1.52 


0.053(55) 


0.0111(6) 


0.15(15) 


0.031(2) 


18.3 


32 


1.52 


0.036(40) 


0.0062(8) 


0.16(18) 


0.028(4) 



Table 6: Results for the damping rates. 

the same with Ca — > 1/2. Let us now assume Co 2 to be this divergent contribu- 
tion, and write Co 2 = 5g 2 T/a. Then S = (2 + 1/2) 0.08606 = 0.2151. Denoting 
the interpolated (3 = 13.6 data by the subscript 1 and the j3 = 21.7 data with- 
out a subscript, we have with ujulg 2 T = 1.14, cluh = 0.210 (cf Table |^) and 
aico H = auo H * (oi/o) = 0.210 * (21.7/13.6), 



UJ\W 



Uljj + UJ 



7,2 



~,2 



Ujfj + Ujf 



+ 5(g 2 T/uj H )/(aujH) 
l + S^T/ou^/ia^H) 



1 + 4.185 
1 + 2.62 5 



1.10, (20) 



for S = 0.2151. This is somewhat smaller than the 'measured' ratio uow/^iw — 
1.61/1.41 = 1.14. Using this as input and solving for 5 would give 5 = 0.38. 
Notice that the 5-terms in the last square root above are not small compared to 
one. In this interpretation all the difference between ujh and low is ascribed to the 
^-contribution, which is probably not the whole story in a more detailed perturba- 
tive calculation, but which works reasonably well. For example, u w /uh — 1 = 0.41 
for the 21.7 data, while ^l + A.18 5 - 1 = 0.38 for 5 = 0.2151. Furthermore, for 
(3 = 12.5, 13, 14, 20 the predicition for u w /u H - 1 is 0.48, 0.36, 0.22, 0.080, 
respectively, to be compared with the data 0.46, 0.36, 0.23, 0.085. The above 
interpretation suggests that the approximate lattice spacing independence of our 
data for the plasmon frequencies is fortuitous and that it may well turn into a 
clear 1/yfa behavior at much smaller lattice spacings. 

6 Autocorrelators in the plasma phase 

In the plasma phase the behavior is qualitatively different. Fig. § shows that 
Cnit) still has roughly the form (|T9"D, but with an average C that gets significantly 
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smaller in a period of oscillation. A reasonable fit to the data could be obtained 
by letting C — * Ct~ s in a region starting a little before the second maximum to 
times where the oscillations appear to have died out, judging by eye. An example 
of such a fit is shown in figure |^. The large statistical errors for Ch in the Higgs 
phase show up in the lack of smoothness in the data and the resulting jackknife 
error on 7^ is very large. The results are in Tables [5], |6] and Figs. ^| and [7[ The 
fitted power 5 is about 0.2 - 0.4. For larger times (t > 100 for 11*) the power 5 
appears to increase, as shown in the double logarithmic plot Fig. |^, which also 
shows that there is no good evidence for exponential behavior. Asymptotically 
one would expect the t~ 3 ^ 2 behavior familiar from zero temperature |36[ . 



The W autocorrelator behaves very differently. Fig. [10] shows what looks like 
Cn(t) in Fig. ||, except that the time scale has shrunk by an order of magnitude. 
The periods of the oscillations are of order 1, which suggests cj's of order 2tc 
and also the 7's appear to be very large in lattice units. There is clearly no 
scaling in this time interval since the period of the oscillation hardly changes as 
/3 increases. We recall that the oscillations start already in the Higgs phase quite 
far from the phase transition, as mentioned above after introducing Fig. |3|. They 
are not analogous to those of C#(£) in Fig. ||, but of Cn(t) in the insert in Fig. |[ 
The small time oscillations correspond to energies of the order of the cutoff and 
must be considered as dominated by lattice artefacts. 

Fig. [11] shows the larger time region on a logarithmic scale. We see good evi- 
dence for exponential behavior. A power behavior appears to be excluded in the 
time range shown (cf. insert). The drop of CV(£) can be fitted to an exponential 
with a damping rate in physical units that is lattice spacing independent, within 
errors, as shown in Table [£] and Fig. [| On the other hand, the slopes of (3 — IV 
and 11* are equal within errors which suggests that the W dynamics is practi- 
cally independent of the scalar field dynamics in the plasma phase. A similar 
effect has been noticed for the screening lengths, see e.g. [p5| . There appear to be 
no oscillations in the exponential regime. At very large times we found irregular 
oscillations with a period of order 600 (18.3 data), indicating a frequency of order 
0.01, but the statistics is insufficient to be able to say anything more definite. 

Ignoring this hint of oscillations at large times as being noise, it is interesting 
to explore an interpretation in terms of overdamping such as discussed in fl6| , [18|j . 
Suppose the W mode satisfies an equation of the form (d 2 + Td t + M 2 )Cw(f) = 
for t > 0, with some basic frequency M and damping rate I\ Substituting 
exp(— iut) gives the equation — uj 2 — iVu + M 2 = 0. For small T 2 <C M 2 this leads 
to damped oscillations, but for large T 2 3> M 2 the system is overdamped, in the 
sense that the solutions uj are purely imaginary and no oscillations occur: u\ ~ 
— iT, u>2 ~ —iM 2 /T. Identifying the smaller \uj2\ with our measured damping rate 
7iy, we can try to interpret Cw(t) for J3 = 18.3 with aT = (9(0.5) (qualitatively 
representing the strong damping in Fig. [Tl], while ignoring the oscillations) and 
aM = ay/jwT = a/0.0062 aT = 0(0.06). However, since the above equation 
for Cw(t) is supposed to hold under the condition \u\ 2 <C M 2 , it seems best to 
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discard u>i altogether. It is interesting that jw is inversely proportional to T. 

7 Discussion 

We have learned that numerical computations of real time correlation functions 
in hot SU(2)-Higgs theory are feasible in the classical approximation. In this 
exploratory study we have not indicated statistical errors in our plots of the 
correlation functions. In most cases these are quite reasonable, but in some they 
would swamp the data lines, e.g. for H in the plasma phase. However, we believe 
the jackknife errors given in Tables [5], || and Figs. ^ [7| are mostly realistic and they 
can be improved with a reasonable increase of computational effort. The errors 
on 7# appear to be overestimated, but this can only be clarified by increasing 
statistics. Systematic errors are harder to assess. It is one thing to compute an 
autocorrelation function, and quite another to extract plasmon parameters from 
these when the damping rate is relatively large. But then the very relevance of 
such concepts is questionable in case of strong damping. 

The results for the H, W plasmon frequencies and damping rates appear to 
be roughly independent of the lattice spacing a. We found, however, that a small 
effect in the W frequency at T/T c = 0.78 could be interpreted as compatible with 
the expected 1/ ' \fa divergence |J"8"| . The approximate a-independence of the H 
frequency may due to the fact that Higgs screening mass is already tuned by the 
counterterms based on dimensional reduction (for the W there is no room for 
such tuning within the classical form of the hamiltonian) . 

As argued in the Introduction, comparison with analytical results in the quan- 
tum theory is to be done in the Higgs phase. The numerical damping rates in 
the Higgs phase compare reasonably well: the (3 = 21.7 (T/T c = 0.78) data 
give j w /g 2 T = 0.22(4) and -f H /g 2 T = 0.028(4), to be compared with 0.176 and 
~ 0.02, respectively [|37|, [38|1 . Of course, these values will depend on the distance 
from the transition at T c (cf. Fig. |G[). The approximate lattice spacing indepen- 
dence of the damping rates may express an insensitivity to the magnitude of the 
plasmon frequency. This is the case in the quantum theory where the plasmon 
frequency is of order gT and the damping rate of order g 2 T, which is the primary 
scale of the classical theory. We note also that the perturbative values for the 
quantum and classical plasmon frequencies for the gauge field are in fact not very 
different in our simulation: for SU(2) gauge theory we have acu^ = ^V^/ft and 



0.8297/ V/? PL which are equal for (3 = 11.6. 



The damping rate associated with H increases substantially when the tem- 
perature is raised above T c , while the rate associated with W drops dramatically. 
Both damping rates appear to be roughly temperature independent in the plasma 
phase, although the errors on 'jn are much too large for a definite conclusion. The 
data for the W correlation function in the plasma phase suggest overdamping, 
with a rather small damping rate. For a proper interpretation we may have to do 
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better than just a pole approximation to the correlation function. It is interesting 
that a similar phenomenon was found recently in a simulation of correlators of 
SU(2) gauge fields, transformed to Coulomb gauge |l||. We stress again that in 



the plasma phase the 'plasmon' frequencies and damping rates of the composite 
H and W fields should not be identified with those of the elementary Higgs and 
gauge fields. 
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Figure 1: Screening masses itih,w / Q 2 T versus T/T c . Solid symbols correspond to 
H, open symbols to W; circle: 24 3 data; diamond: 20 3 data; square: 32 3 data. 
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Figure 2: Scan of C#(£) from the 24 3 data. The first minima in t > 5 are in order 
of decreasing J3 = 20, 14, 13, 12.5, 12. The lines in the middle of the Ch plot are 
for (3 < 12, shown enlarged in the insert. 



19 



0.04 - 



0.02 







■0.02 



"i i i r 



"i r 



"i r 




j L 







5 



10 
time 



15 



Figure 3: Scan of Cw(t) from the 24 3 data. The first minima in t > 5 are in 
order of decreasing /3 = 20, 14, 13, 12.5, 12. The lines ending before t — 10 are 
for (3 < 12; bottom to top (t — 4): (3 — 11 (touching the minimum of (3 — 12), 
10 (nearly coinciding with 11), 8, 6 and 4. 
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Figure 4: The autocorrelator Cn(t) in the Higgs phase for (3 = 13' (small ampli- 
tude) and 21.7 (large amplitude). 
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Figure 5: The autocorrelator Cw{t) in the Higgs phase for j3 = 13' (lower tail) 
and 21.7 (upper tail). 
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Figure 6: Plasmon frequencies w/g 2 T versus T/T c . Solid symbols correspond to 
H, open symbols to W; circle: 24 3 data; diamond: 20 3 data; square: 32 3 data. 
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Figure 7: Damping rates , y/g 2 T versus T/T c . Solid symbols correspond to H, 
open symbols to W; circle: 24 3 data; diamond: 20 3 data; square: 32 3 data. 
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Figure 8: Left: Cn(t) in the plasma phase; /3 = 18.3, 11* and 11' in ascending 
order. Right: fit (dashed curve) to the data at (3 — 18.3; the fitting region 
corresponds to the data shown. 
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Figure 9: InCnit) versus t (left) and lnC#(£) versus lnt (right) in the plasma 
phase. 
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Figure 10: Short time behavior of Cy/{t) in the plasma phase. Top to bottom 
(t» 1): (3= 11', 11* and 18.3. 
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Figure 11: Plot of lnCV(£) in the plasma phase for larger times; (3 = 11' (top), 
11* (middle) and 18.3 (bottom). The insert shows the same time interval on a 
double logarithmic plot. 
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